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^ ! Abstract 

' We examine the relation between Coulomb-gauge fields and the gauge-invariant 

fields constructed in the temporal gauge for two-color QCD by comparing a vari- 
ety of properties, including their equal-time commutation rules and those of their 
conjugate chromoelectric fields. We also express the temporal-gauge Hamiltonian 

■ in terms of gauge-invariant fields and show that it can be interpreted as a sum of 
the Coulomb-gauge Hamiltonian and another part that is important for determin- 

I ing the equations of motion of temporal-gauge fields, but that can never affect the 

■ time evolution of "physical" state vectors. We also discuss multiplicities of gauge- 
invariant temporal-gauge fields that belong to different topological sectors and that. 



o , 

O , in previous work, were shown to be based on the same underlying gauge-dependent 

O ' temporal-gauge fields. We argue that these multiplicities of gauge- invariant fields 

are manifestations of the Gribov ambiguity. We show that the differential equation 
, that bases the multiplicities of gauge-invariant fields on their underlying gauge- 

^j' dependent temporal-gauge fields has nonlinearities identical to those of the "Gribov" 

^ ■ equation, which demonstrates the non-uniqueness of Coulomb-gauge fields. These 

multiplicities of gauge-invariant fields — and, hence, Gribov copies — appear in the 
k>( \ temporal gauge, but only with the imposition of Gauss's law and the implementa- 

^ ' tion of gauge invariance; they do not arise when the theory is represented in terms 



of gauge-dependent fields and Gauss's law is left unimplemented. 



1 Introduction 

In earher work, we have implemented the non-Abelian Gauss's law that applies in QCD 
by constructing states that are annihilated by the "Gauss's law operator" ^'^(r) for the 
temporal {A^ = 0) gauge, |l| where, for two-color QCD, 

^"(r) = d,U'^{r)+ge^'^A%T)U'^^ir), g%r) = g%r)+f,{r), and f,{r) = g ^\r)^i^{r) 

(1) 

and where n"(r) is the negative chromoelectric field as well as the momentum conjugate 
to the gauge field Aj{r). We have, furthermore, used the gauge- invariant quark and gluon 
operator-valued fields constructed in Ref. to transform the QCD Hamiltonian into a 
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form in which it is expressed in terms of these gauge-invariant fields. |^ Most recently, 
we have studied the relation of the gauge-invariant to the gauge-dependent gauge fields 
in the temporal gauge. In particular, we have solved the nonlinear integral equation 
that expresses the requirement that the non-Abelian Gauss's law be implemented, and 
have discussed the consequences that the solutions of this integral equation have for the 
topology of the gauge-invariant gauge field. [|] 

In this paper we will address a number of questions that pertain to the relation between 
QCD in the Coulomb gauge and our formulation in which QCD in the temporal (Weyl) 
gauge is expressed entirely in terms of gauge-invariant fields — i.e. in terms of operator- 
valued fields that commute with the generator displayed in Eq. (0). QCD in the Coulomb 
gauge has been discussed by Schwinger, ||^, ^ by Christ and Lee, 0, |^ by Creutz et. al., 
and by Sakita and Gervais, []TU| among others. Quantization of QCD in the Coulomb gauge 



encounters a number of difficulties: These include the well-known Gribov ambiguity that 
becomes an impediment when Gauss's law is inverted to solve for the timelike component 
of the gauge field. [|ll|, |12| Attention has also been called to operator ordering problems 
encountered when noncommuting fields appear multiplicatively in the Coulomb-gauge 
QCD Hamiltonian. 0, ^] Some authors have circumvented the latter problem by treating 



the Aq = gauge fields as a set of Cartesian coordinates and the Coulomb-gauge fields as 
a set of curvilinear coordinates and using standard methods to transform from the former 
to the latter. |] 

We will show in this work that the gauge- invariant fields we have constructed — the 
gauge fields as well as the chromoelectric field — have properties that so closely match 
those of the corresponding Coulomb-gauge fields, that they can be identified with them. 
The gauge-invariant gauge field we have constructed is transverse, as is the gauge field in 
the Coulomb gauge. The gauge-invariant fields obey commutation relations that are the 
same as those in the Coulomb gauge with the exception of operator order in the commu- 
tator. We will, in this work, construct a number of operators that are important in QCD 
dynamics — the Hamiltonian, the Gauss's law operator, the gluon and quark color-charge 
densities, etc. — so that they are expressed entirely in terms of these gauge-invariant 
fields (for short, we will call them the "gauge-invariant" operators). We will show that 
this gauge-invariant Hamiltonian is not the same as the Coulomb-gauge Hamiltonian, but 
that it contains a version of the latter accompanied by another term that is specific to the 
temporal gauge. This additional term in the gauge-invariant temporal-gauge Hamiltonian 
has no dynamical consequences in the subspace to which the physical state variables are 
restricted, so that the two operators — the Coulomb-gauge Hamiltonian and the gauge- 
invariant Hamiltonian — have identical physical consequences. However, the additional 
term in the gauge-invariant temporal-gauge Hamiltonian affects the equations of motion 
of the gauge and quark fields. Finally, we will also discuss the relation between the mul- 
tiple solutions of the equations for the gauge-invariant fields that we obtained in Ref.Q] 
and the Gribov copies of the Coulomb-gauge fields. 
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The plan for this paper is as follows: In Section we will compare the commutation 
rules for the Coulomb-gauge field with those for the gauge-invariant field that we con- 
structed in the temporal gauge, and discuss the sense in which these two fields can be 
identified with each other. We will also make a similar comparison for the gauge-invariant 
momentum (and negative chromoelectric field). In this section, we will also construct a 
number of gauge-invariant operator-valued quantities, such as the gauge-invariant gluon 
and quark color-charge and color-current densities, and the gauge-invariant version of the 
Gauss's law operator given in Eq. (|^). We will, furthermore, express the temporal-gauge 
Hamiltonian in terms of these gauge-invariant variables, and discuss its relation to the 
Coulomb-gauge Hamiltonian. In Section ^, we will discuss the multiple solutions of the 
nonlinear equation that determines the gauge-invariant field, and their relation to the Gri- 
bov copies of the Coulomb-gauge field. In Section ^ we will conclude with observations 
based on results presented in earlier sections. 



2 Gauge- invariant temporal-gauge fields 

The quantization of QCD in the temporal gauge avoids many of the problems encountered 
in quantizing QCD in the Coulomb gauge, but at the expense of leaving Gauss's law still 
to be implemented after the quantization has been carried out. In quantizing QCD, we 
make use of the Lagrangian for two-color QCD, in which the gauge fields are in the adjoint 
representation of SU(2):[] 

L = -\f^F^ + + Jt^ - Jo^S + ^gauge - ^P {m - ■ d) ^j; (2) 

where 

Ft^ = d,Nl - d^A] - ge-''^A\A] , F^^ = doA^ + + ge'^'"^AtA^, , (3) 

with jf = #"faiy^, and j^ = g^^—^- (4) 

-^gauge is a gauge- fixing term. When the gauge-fixing term —AqG"" is used in Eq. (|), 
and the Dirac-Bergmann method of constrained quantization is used, [O, IT1| the Lagrange 



multiplier field is incorporated into the time-derivative of IIq, which, in this case is 
Djll" -|- j'q + G"^. The presence of the Lagrange multiplier field in the secondary constraint 
DiU'^+j^+G"' = terminates the chain of secondary constraints very quickly, and leads to 
Dirac commutators that differ in only trivial ways from canonical Poisson commutators. 
However, the Gauss's law constraint is not imposed in that process. Alternatively, it is 
possible to entirely avoid the need to consider primary constraints in the temporal gauge 

"'^we use nonrelativistic notation, in which aU space-time indices are subscripted and designate con- 
travariant components of contravariant quantities such as Af or jf, and covariant components of covariant 
quantities such as di. Repeated indices are summed from 1^3. 
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by using the gauge-fixing term — SoAgC^ instead of —A^C^, so tliat —C^ becomes tlie 
momentum canonically conjugate to Aq. The gauge constraint then is OqAq = 0, which, 
with the imposition of = and Gauss's law at one particular time, implements both. 



the gauge condition and Gauss's law for all times; 1 15] the same procedure can be extended 



to all axial gauges for which the gauge condition is y4Q + 7y43 = 0, where 7 is a variable real 



parameter. |T6| Even in these cases, however, the implementation of Gauss's law at one 
time — at t = 0, for example — is still necessary. Finally, a very direct way of quantizing 
QCD in the temporal gauge is simply to set = in the original Lagrangian. Gauss's 
law is not one of the Euler-Lagrange equations in this formulation, and must be imposed 
after the basic quantization has been carried out.|0, IT8|, [19 



Ref. [Q addresses the imposition of Gauss's law in the temporal gauge by explicitly 
constructing the states that are annihilated by the non-Abelian Gauss's law operator 
given in Eq. ([^). The mathematical apparatus required for that purpose also enables 
us to construct the gauge-invariant gauge and quark fields. This apparatus includes the 
defining equation for a nonlocal operator-valued functional of the gauge field — the so- 
called "resolvent field" A]{r) — which has a central role in this construction. In the 
two-color SU(2) version of QCD, with which we are concerned in this work, the resolvent 
field appears in the gauge-invariant gluon field in the form 



[^Gu(r) ^ ] = Vc{r) [A'^{r)-4]Vc\r) + ^^Vc{r) d.Vc\r), (5) 
where Vc{r) incorporates the resolvent field, as shown by 

Vc{r) = exp ( -z(?3^(r)f ) exp (-z(;A'°(r)f ) (6) 

and ^^^(r) = exp (z(^A'"(r)f ) exp [tgy^{r)^) (7) 

with X'^ir) = §Af{r) and 3^(r) = |:4f(r). (8) 
The composition law for two successive rotations can be used to express Vc{r) in the form 

K:(r) = exp(-z(72-(r)f) (9) 



where Z°'{r) is a well-known functional of A:'"(r) and 3^"(r). The transversality of the 
gauge-invariant field is manifested most directly by transforming Eq. into 

4u(r) = A^.(r) + [5., - f^]:4f(r), (10) 

demonstrating that Ag|j(r) is the sum of A^^{r) — the transverse part of the gauge field 
A^(r) in the temporal gauge — and the transverse part of the resolvent field. The resolvent 
field is also required for defining the gauge-invariant quark field 

V^Gi(r) = Vc(r)^(r) and ^^^(r) = ^t(r) i/-i(r) . (11) 
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Vc{r) can be understood as an extension, to non-Abelian gauge theories, of a much simpler 
but similar operator that made charged QED states gauge-invariant and that was intro- 
duced by Dirac.|2^ Dirac-like operators implement gauge invariance without introducing 
path dependence, and it has been argued that such Dirac-like operators have advantages 
for implementing gauge invariance over other procedures that generate gauge-invariant 
charged fields with path-dependent line integrals. ^ We have noted that Eqs. 



and (0) implement gauge invariance — they do not describe gauge transformations. fT]] 
By construction, Vc(r) transforms A" and ip so that they become invariant to further 
time-independent gauge-transformations consistent with the temporal gauge condition.^ 
Although formally Eqs. (||) and (|Tl|) are very suggestive of gauge transformations of gauge 
and spinor fields respectively by c- number functions, that is not what they are. Vc{r) is 
itself a functional of gauge fields, and is transformed by the same gauge transformations 
that transform the fields on which it acts, just as is the case in the corresponding operator 
Dirac proposed for QED. We will show, in this section, that Eq. (^) describes a gauge 
field that can be identified with the Coulomb-gauge field; that, in combination with other 
suitably constructed gauge-invariant quantities, it can be used in the representation of the 
temporal-gauge Hamiltonian; and that the dynamically effective part of this Hamiltonian 
can be interpreted as the Coulomb-gauge Hamiltonian. 

Eqs. (^-(|n|) show that the resolvent field ^^(r) has a central role in the representation 
of QCD in terms of gauge- invariant field variables, ^^(r) is determined by a nonlinear 
integral equation that was obtained in the course of constructing the states that implement 
the non-Abelian Gauss's law.[|I| In subsequent work, this nonlinear integral equation was 
transformed (subject to an ansatz) into a nonlinear differential equation that we solved, 
resulting in nonperturbative representations of the resolvent field and the gauge-invariant 
gauge field. ^ We will discuss these nonperturbative solutions in more detail in section |^. 
An important corollary of the formalism that leads to this equation for the resolvent field 
is the commutation rule |2B 



J dr [u%y),Alir)] W,T(r-x) + (y-x) = 



2 r=0 



where W,7(y - x) = 5,„Wfe,(y - x) (13) 
with Wfe,(y-x) = -z('5,fe-^')5(y-x). (14) 



^Gauge transformations within the temporal gauge, for which the Gauss's law operator given in Eq. 
(P is the generator, are restricted to time-independent gauge functions, so that the temporal-gauge 
constraint is not violated. 



5 



In Eq. (|T^, we have used a notation introduced in earlier work: ^, ^] ef'"' represents 
the chain of structure constants 



^vba 



gVlbsi ^SlV2S2 ^S2VZS3 



g«{r-2)'"(r-l)«(r-l)g'5(r-l)''ra 



(15) 



where repeated superscripted indices are summed from 1^3. For r = 1, the chain reduces 
vha _ ^vba. „„A f^v ^ — n ^vha _ _^ Furthermore, we have used 



to 



e"™; and for r = 0, 



(16) 

which, for r = 1, reduces to ^ (7^f)^(x)Wfci(x)) = || (A^, fc(x)Wfci(x) ) and for r = 0, 
reduces to ^ (^T^q) ^(x)Wfcj(x)^ = Wjj(x). For clarity, we give the explicit form of the r-th 
element of J {T^^^ f^{y)Uki{y - x)), which is: 

(iz(2) 



92 



{^)kiJ)UM{y-^))=^{-lr^j 



dyjJ 47r|y - z(i)| 



9^z,(i)y 47r|z(i)-z(2)| 



d 



dzi, (2) 



(izfr-l) 



-AJi^ fx + Mz 



47r|z(r-2) — z(r-l' 
1 



d 



V47r|z(r-i) - x| ^'^ 
the leading (0-th) order term is 



4'- 



47r|z(r-i) - zf'^"=^^''9z^9z^47r|z-x| 



X 



-2 



92 



(17) 



(1^ 



We observe that U^Ky — x) serves as a projection operator that selects the transverse 
parts of the fields over which it is integrated, and that Eqs. and ([l2|) in combination 
show that 



dr 



n?(y),^^,(r) W-(r-x)+W?r(y-x)= n;(y),A^„(x 



(19) 



The trace Tr[rVc"V"\/c], which appears on the right-hand-side of Eq. (|12D, can be ex- 
pressed as {T'')qn{Vc''^T'^Vc)nq] with the use of the identity 



b b 
' kI qn 



we obtain 



Vc-'t'Vc 



{r'')AT'),r,{Vc-'T'Vc)nq = 2 (K:"'^'^:)^^ - Tr 
Since Jr[Vc'^T'Vc] = Tr[r"] = 0, it follows from Eqs. ([l|) and (||) that 

'_{T')^U%y),AlM] = Y: {Vc\y)T^Vc{y)) (_i)'-+V4-| (7^f),(y)W,7(y-x 



(20) 
(21) 



r=0 



(22) 
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and that, upon contraction with [Vc(y)T'^V^ '^{y)]iK, we obtain 



n^,,(y),A^,,(x) 



r=0 



djdk 

Q2 



5(y - x) 



where 



(23) 



(24) 



has been identified in previous work as the gauge-invariant momentum conjugate to the 
gauge field, (and the negative gauge-invariant chromoelectric field). P] The following fur- 



ther explanatory remark for Eq. 



can be given: We let IIgi i = Hq 



GU~5 ^Gl 



^GU 2 ' 



Ai = ^, and make use of the analogy of Eq. (|g) to a gauge transformation to define a 
corresponding Aq\ q, and to observe that, since in the temporal gauge Aq = 0, 

Aoi = ~^Wo where Wq = VcdoVc^ . 

We recall an operator-order we have previously found necessary, and impose it on the 
obvious definition of a "gauge-invariant momentum (and negative chromoelectric field)" 



IIgu = WdiAoi + doAoii - ig [Agm, Agi o] 



(25) 



where, using a notation introduced in Ref. bracketing between double bars denotes a 
normal order in which all gauge fields and functionals of gauge fields appear to the left 
of all momenta conjugate to gauge fields, but where that order is imposed only after all 
indicated commutators have been evaluated (including the commutator implied by the 
derivatives do and di). We observe that 



I^O^GIi 



VcA,Vc-\ Wo] + VedoA^Vc' + {doVcdiVc' + Vcdod.Vc') 



and 



I|MgioII = -| 

- ig [^Glh^Glo] 



^^Vc^oVc~' + VcdodiVc-' 



VcAiVc\W, 



Combining Eqs. (E6|-B3), we find that 



n 



Gh 



\Vrd.A.Vr-^\ 



doVe.diVc 



Vc^Vc-'Ul 



-1 



(26) 

(27) 
(28) 

(29) 



which agrees with Eq. (0). This demonstration is not necessary to prove that is 
gauge-invariant — the transformation properties of 11" and of V^yV^"^ suffice for that 
— but it nevertheless makes the identification of IIq ■ as the gauge-invariant momentum 
conjugate to AQ^■ more understandable. 
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As we will point out later in this section, the commutator given in Eq. ( PB]) agrees 
with the corresponding commutator for the gauge field and its conjugate momentum in 
the Coulomb gauge, given by Schwinger. The gauge field in the Coulomb gauge and 
the gauge- invariant gauge field AQ^^{r) constructed within the temporal gauge are both 
transverse.0 In Ref. |^ and in our work, the momentum conjugate to the gauge field — 
the negative chromoelectric field — has a longitudinal component, which is needed to 
implement Gauss's law. The Coulomb-gauge commutation rules in Refs. H, |^ differ 
from Eq. (pSj), because in these works, the momentum conjugate to the Coulomb-gauge 
field has been defined to be transverse. 

We next turn to a discussion of the "gauge-invariant Gauss's law operator" — the 
Gauss's law operator in which gauge-invariant fields replace the original gauge-dependent 
ones — which we define as 

^G,(r) = d.Ui,{r)+ge'^^Al,,{v)Ul,,ir)+f,M ■ (30) 

In particular, we want to investigate whether the use of the gauge-invariant Gauss's law 
operator as the generator of gauge transformations is consistent with the gauge invariance 
of y4Qj(r). We will demonstrate this consistency in this section. 

We will first demonstrate a simple relation between ^gi(^) ^'^(i")) from which 

dr[gl{r),Al,,{^)] 6uj\r) = (31) 

is an immediate consequence. We observe that 

d,Ui,{r) = lT,{r%Vc{r)r%-\r)}u'i{r) + '^Tr {r^%{r)r%Vc-\r)}u^{r) + 

iTr{rVc(r)rVc-i(r)}9A'(r) (32) 

and that, for Xi(r) = Vc{r)diV(r^{r), 

d.Ui, = iJr {r^ [Vcr%-' , x^] } + ^Tr {t'^Vct'Vc'} d^U^ (33) 

In an Appendix, we will show that we can set 

Xi = t''^ui and Vct'Vc-' = T^'n^t, (34) 

where Vui and TZ^b are functions of gauge fields only — they are independent of the 
canonical momentum 11^ and also contain no further SU(2) generators — so that 

iJr {r^ [Vcr'Vc' , x^] } = fTr {r' [r^ , r"] } 7^.,P„.^^ 

= e'^™7^,bP„,^,^ (35) 



•^The relation between gauge-invariance and transversality can be inverted to produce a perturbative 
representation of the gauge-invariant field, as in Refs.[pT|, p^ . 
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and 



Similarly, using Eqs. (^), (p^), and (p4D, we can express ge'^'^^ Aq\^Wq\^ as 



(36) 



(37) 



and find that e'^'^'"lZvbVui^\ terms in Eqs. (^) and 
Gauss's law operator can be expressed as 



cancel, so that the gauge- invariant 



where, in order to obtain the last term in Eq. (pH]), we have usedQ 

JoGi = g'^^Vc^r^Vc'ilJ = -JZdhgi^W'i: = 7^,,Jo' 
which we justify in the Appendix. Eq. (0) leads to 

7^d6 = '^Tr[T''Vcr%-'], 

so that 

e'^"^7^„,7^,, = le'^-^'ir^Ur^kp (Vct'^Vc-') .. (Vcr%-' 
We make use of the identity 



pk 



e )kp 



to obtain 



[Sjk{T'^)ip - 6ip{T' )kj 



T ,r 



and after relabeling dummy indices, observe that Eqs. ( p8| ) and ( ^3]) lead to 



(38) 

(39) 

(40) 
(41) 

(42) 
(43) 

(44) 



Since TZab trivially commutes with AQ^^, the fact that commutes with is sufficient 
for the demonstration that also commutes with Aq,, thus proving Eq. (|3TD. 

Another, very direct demonstration of Eq. (|3T| ) begins with the use of Eq. (^3]) to 
obtain 



dr 



9,tfG„(r),^G,.(x) 



r=0 



(45) 



*To simplify this argument, we make use in this discussion of a representation defined in Ref. Q — 
the so-called C-representation — in which Vc(j)ip(r) is the gauge-invariant spinor and in which Jq^, = 
5^^V^~^(t°/2)VcV' is the gauge-invariant color charge density. Elsewhere in this paper, we have used the 
so-called A/'-representation, in which the Vc transformation has already been implicitly carried out for 
the quark field (but not the gauge field), so that it is the quark field ■0 and Jq = g-ip^ {t°' /2)^ that are 
gauge invariant. 
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we observe that the r = term on the right-hand-side of Eq. ( ^31 ) vanishes, because the 
degenerate values hsted immediately preceding and following Eq. ([T6|) demonstrate that 
9j(^7^Q-) ^(r)Wfcj) = 0. We can therefore begin the sum with r = 1 instead of with r = 0, 
redefine the dummy index r to be r + 1, and then initiate the sum with r = for the new 
index r, obtaining 



(46) 



r=0 



We can also evaluate 



gidye'^-Al^^iy) ln^,^.(y) , A^^uW] 5u\y) = I dyY^i-iy^'g^^'e''^^^^'^^ 



r=Q 



^G„(y)^ (^^^).(y)w..(y-x)) Suj\y) 



and observe that 



d 



e'^-erAl,^{y)^ {%,{y)Uu.{y-^)) = e^, (7^f+,),(y)W,,(y-x; 



(47) 



(48) 



so that 



gjdye'^-Al,^{y) [n^,,(y) /1"gu(x)] 5uj\y)=j dyY.{-ir'g'^'4% (^f+i) .(y)W..(y-x)) 5u\y). 

r=0 

(49) 

Cancellation between the right-hand sides of Eqs. (^) and (|^) verifies Eq. (^l]), and 
therefore confirms that the use of as the generator of infinitesimal gauge transforma- 
tions is consistent with the gauge invariance of j(x) (but not with the gauge invariance 

ofn^u(y)). 

The mathematical apparatus we developed for constructing gauge-invariant fields en- 
ables us to express the QCD Hamiltonian entirely in terms of these gauge-invariant fields 
and n^i^. (the gauge-invariant negative chromoelectric field). [0, ^ The QCD Hamiltonian, 
represented in this way, has the form 



Gl 



Hq\ + Hq 



(50) 



where Hg annihilates states that implement Gauss's law and has no dynamical conse- 
quences in QCD. Hq\ is the effective Hamiltonian in this representation of QCD and is 
given by 



Gl 



dv 



|H"G,Ur)n^,,(r) + \F^,,,{v)F^,,,{r) + V^t(r) (/3m - laA) V^(r)] + H' , (51) 

where, in this case, ip and ip'^ denote the gauge-invariant spinor (quark) fields, and where 
F^,,^{v) = d,Al,,{r) - d,Al,^{r) - ^7e'^^^A^„(r)A^,^.(r). (52) 
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H' is given by 

H' = jdv {\j;U{v)^,lCl{v) + i/Co"(r)ljo%,) - i/CS(r)l/Co"(r) - jT(r)A^Gu(r)) 

(53) 



where Jo%,) = ^^e'^^'^A^GuWn^GuW (54) 

is the gauge- invariant glue color-charge density, jf(r) = gilj\v)ai^ip{Y) is the gauge- 
invariant quark color-current density in this representation, and where 

/Co'(r) = E 4n) (-1)"^?" (^J)(r)jo'(r)) , (55) 

With Tl,{v)m = 4« ,)(r) %i «) ,)(r) ^ (■ ■ ■ (^^^^(r) %I (j-(r))))) . (56) 
The more transparent explicit form of /Co(r) is: 

/C^(r) = -Jo^(r) +,e^--<-(r)|- / ^^^Jo(x) + 

d f rfy(i) ^ 



47r|r - y(i)| 



n—2) 4 t'(n-i) / N 



dy{n^3)iJ 47r|y(„_3) -y(„_2)| ^ ^l/(n-2) i 



/ 



-^5^^^=^^ r^cS (y..-.,) / .., ;S(x) + ■ ■ ■ . (57) 



47r|y(„_2) - y(n-l)| <9?/(n-l) A: 47r|y(„_i) - X 



In formulating Eqs. (|5TD -(^7D, we note that while it is trivial that AQ^■{r) is hermitian, and 
while we will show that QQ\{r) is hermitian as well, IIq -(r) and IlQ^■{r) are not identical. 
As can be seen from Eq. (|2^), in order for ng|^(r) to be hermitian, n^(r) would have 
to commute with -R(if,(r) (defined in Eq. (ppj)). We therefore distinguish between IIqi j(r) 
and nQj(r) in this discussion; similarly, Eq. (|23| ) shows that Jq and Jg ^ will differ. The 
Hamiltonian and its components, Hq\, Hq\, Hg, and H', are all manifestly hermitian. 

Eq. (^3|) is very suggestive of an interaction Hamiltonian in the Coulomb gauge, but 
with /Cq appearing in the position in which one might expect to find the charge density 
of the matter field j^. )Cq contains Jq as a crucial component, but extends the latter 's 
dynamical effect over a greater region in space than j'q itself occupies, through a series of 

similar point is made in Ref. note, however, that the chromoelectric field in Ref. |^ is transverse, 
while our chromoelectric field includes its longitudinal part, as in Ref. the longitudinal part of the 
chromoelectric field accounts for its lack of hcrmiticity. 
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"chains" of interactions in which each hnk has the form e"''"^ Aq^^. The effect of these 
chains of interactions is that if, for example, j'q describes quark color-charges that are 
limited to a relatively small volume, /Cg could be significant over a substantially larger 
region of space. Further investigation of these nonlocal interactions can be facilitated by 
the observation that 

IC^n + Qe''''''Aa~n = -j^. (58) 



Q2 



The issue here is not that JCq is gauge-invariant. In fact, in the representation used to 
express H', both /Cq and Jq cire gauge- invariant, as has been discussed in Ref. 

A further significant feature of H' is that, besides the nonlocal interaction — |/Cq^/Cq, 
it also includes additional interactions of the gauge-invariant glue color-charge density 
Jq with /Cq, so that quark and glue color charge-densities interact with each other in 
the non-Abelian theory in which the gauge field, as well as the matter fields to which it 
couples, carry color-charge. 



Hg — the part of the Hamiltonian that annihilates states that implement Gauss's law 
and therefore has no dynamical consequences — can be given as 



He 



-\ I dr 



^a^^oW + /Co"(r)^^a, 



(59) 



Since is hermitian, and since any state |\E') for which Qq\{^) |\1') = will time-evolve 
so that ^q(x) exp(— ii^Gi ^)|^) = as well — properties of that we will prove later 
in this section — Hg will not contribute to matrix elements for physical processes or 
have any affect on the properties of states in the space of "physical" states which must 
implement the non-Abelian Gauss's law. 



To demonstrate the hermiticity of G^ir), we use Eq. (|4|) to show that 



A.(r) 



(60) 



From the proof that Vc(r) transforms gauge- dependent temporal-gauge quark and gauge 
fields into gauge- invariant fields as shown in Eqs. (|) and (pJ|),p5| we observe that 

my) , ^c(x)] = (7Vc(x)^5(x-y) and \g^{y) , ^^"^(x)] = -<7^17^(x)5(x-y) , (61) 



so that 



(62) 



A,(r) = iTr {r' [g%r) , Vc{r)T^Vc\r)] } = 

and ^Q(r) is shown to be hermitian.0 

^We interpret the commutators [C/'=(r) , Vc(r)] and [^^(r) , V^~"^(r)] as limr'_r[^'^(r') : ^(i")] ^^^^ 
limr'^r[S'^(r') , V^~^(r)] respectively to regularize the argument of the delta-function. 
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We will examine the hermiticity of nQ^(x) by first evaluating [nQ|^(x) ,ng|.(y)] and 
comparing it with the commutator [IIq -(x) , nQ^(y)]. In this way, we avoid the necessity 
of evaluating the singular commutator [n^(x) ,i?fea(x)]. To evaluate these commutators, 
we first replace jo{r) by —g^6{r — x) in Eqs. (7) and (19) in Ref.pf, a step justified by 
the fact that both sets of quantities obey the same closed commutator algebra. In this 
way, we obtain 

yc(x)n5(y)v,-i(x) = n5(y) + y: g^'i-iTRMi''^ (^l)(y)^?4'5(y - x)) (63) 

n=0 

and therefore that 

n;(y) , Vc{^)r'^Vc\^)] = ^ ^7"+^(-l)"+^i?..(y)eM {^!)(y)5(y - x) , Vc{^)t^Vc\^)] } 

n=0 

(64) 

from which 

n;(y) , RU^)] = I Y: g^^\-lT^'RMl'^e^''^% [rUy)5{y - x)i?,„(x)} (65) 

follows; and, using Rcb{j)Rdb{y) = ^cd, we also obtain 

ng,,(y) , RU^)] = ^ E ^7"+^(-i)"+^ere'^^"| {rUy)S{y - x)i?..(x)} . (66) 



n=0 



Since 



n^„(x),ng,^.(y)] = [n^,,(x),i?^,(y)]n5(y)- [ng,^.(y),i?„,(x 
it follows that 

naH(x),n^„(y)] =^E^?"^'(-lr' 



n=0 



ef'^e^^^f {^'„)(x)5(x-y)n^,,.(y)} 
-er6'^^"|{7^1)(y)5(y-x)n^„(x 
The explicit form of the n-th order term ^ |7^^)(y)5(y — x)!!^! j(x)| is 



47r|z(i) - z(2)|' 



d 



1 



dzi 



47r|z(n-i)) - z( 



X 



n2u(x) 



dz{n)i^ 47r|z(n) - x| 
and the leading {n = 0) term of the commutator given in Eq. (p7| ) is 

I „^ , . d 1 



a/37 



^a^47r|x-y| ^'^'^^^ ^ %47r|x-y|''^' 



nz,. X 



(67) 



(68) 



(69) 
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The delta-functions that appear in Eq. ehminate the integrations over the last of 

the inverse laplacian in the chain described in Eq. ([561). When there is no delta-function 



in the expressions on which 7^^^ acts, the last inverse laplacian is also integrated over, as 
can be seen by comparing with Eqs. (|55| ) and (^7|). 



We can apply the same procedure used to obtain Eq. (|68| ) to the commutator of 
ng|j(x) and the hermitian adjoint of IIq (y), in which case we get 



x),n^,Vy) 



Raa{^)u%y) [n^x) , RM] - n^„(y) , Raa{^) n^(x) , (70) 



and therefore that 



TTtt 



;x),<Uy) 



n=0 



+1/ 



\n+l 



- ere'^^-l {rl){yny - x)mGu(x)} 

- [RM^^n^)] fi?.a(x),n;(y) 



(71) 



An alternate expression for Eq. (^) can be obtained by defining P(x, y) as an inverse of 
the Faddeev-Popov operatoif] {D-dY^ = {6ab + ge''''''Al^^^)d'^, given by 



Z}-9(x)P(x,y) = 5(x-y). 
15 (x, y) can be expanded as a series in the form 

^^''(X, y) = E 4n) (-1) V7^n)(x)5(x - y) 



(72) 



(73) 



ra=0 



^2 

SO that Eq. (^) can be expressed as 

nau(x) ,n^,^.(y)] = z^? {9,P°'^(x,y)e'^^^n^,/y) - d,V^'{y,^)e'^^Ul,{^)} . (74) 

In the form given in Eq. (|7^), the equal-time commutation relation in Eq. (^) can be 
seen to be in agreement with the one given by Schwinger in Ref. modulo the fact that 
Schwinger's Coulomb-gauge chromoelectric field operators are defined to be hermitian, 
whereas ours are not. The operator-ordering in the commutators in Ref. [^] that corre- 
spond to our Eq. (^) therefore differs from ours. We have defined the gauge-invariant 
chromoelectric fields so that the Hamiltonian in Eq. ( pTD has a simple and tractable form. 
Furthermore, with D^"- = Sbadi + ge'^°'AQ^^, Eq. (f73D, and the identity 



's,,Sba + I E 4) (-l)"^?"^f)(x)A'"(x)) 5(x - y) 



E 

n=0 



^Sbaf 
-'{n){ 



^ik — 



didk 



Q2 



5(x 



(75) 



'''In this, as well as in a number of similar cases, the gauge-invariant field ^ replaces the standard 
gauge-dependent temporal-gauge field in the gauge-covariant derivative. 
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5]) to substitute explicit expressions for the commutators in the 
but comparison of Eqs. (|68D and (|7TD suffices to confirm that 



we can also confirm the agreement between Eq. (^) and the corresponding commutation 
rule given by Schwinger in Ref. 0, again modulo operator order in the corresponding 
expressions. Our commutators in Eqs. (|23| ) and ( |68D differ much more substantially from 
those in Refs.J^, ||, Q because, in these references, the Coulomb-gauge chromoelectric field 
is defined to be transverse, in contrast to the practice followed in Schwinger's treatment 
of the Coulomb gauge as well as in our work. 

We can use Eq. 
last line of Eq. _ _ 

^G\jiy) differs from nQ^(y). It is also possible to evaluate IlQ^j{y) — IlQ^j{y) directly. 
The resulting expression is clearly nonvanishing, but the expression is lengthy as well 
as singular and it is not necessary to quote it in detail in order to make the point that 
the gauge-invariant momentum (and negative chromoelectric field) is not hermitian. We 
observe that ^jllQi • and g^'^'''^AQ^■IlQ^■+jQQ^, the two component parts of ^q, are separately 
gauge-invariant, but are not separately hermitian. 

It is also of interest to examine whether the commutator algebra of the gauge- invariant 
Gauss's law operators closes, by examining the commutator [^^(x), ^Gi(y)]. From Eq. 
we observe that 



and, using 



m^), UMiy)] = ige""'n,,{y)5{^ - y). 



that 



^(?e«^^^a,(x)5(x-y). 



(76) 

(77) 
(78) 



This shows that the commutator algebra of the gauge-invariant Gauss's law operators 
closes and, in fact, is almost identical to the algebra of the gauge dependent Gauss's law 
operators except for a relative sign change between the two cases. 

One consequence of the time-independence of the non-Abelian Gauss's law operator 
given in Eq. (P is the fact that [if, ^°(x)] = 0, so that when a state implements Gauss's 
law (i.e. when ^"(x) |\E') = 0), the time-evolved state ^'^(x) exp(— ii7t)|\&) also vanishes. 
It is therefore important to inquire whether a similar property can be ascribed to the 
gauge-invariant non-Abelian Gauss's law operator ^q, so that the set of states that im- 
plement ^Gi(x) 1^) = also implement that same constraint at all later times, when 
exp(— zi^Gi^) is the time evolution operator. We address this question by observing that 



^a(x),tfG,(y)l =E 



r=0 



'^'^|(^f)(y)5(x-y)) ^^,(x) 



(79) 



follows from Eqs. ( |66D and (^4[) , and that (x) commutes with 1Cq{y) and with Jo(r). The 
commutator of Hq\ and (x) therefore receives contributions only from commutators of 
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nGi(y) and ^gi(x); and, when [n^|(y), ^q(x)] is to the left of another operator, as in 
/rfy[n^Gi(y), G'ai^mUy) or in Jdyge'P'^Al,{y)[UUy), g^^,{x)]{d-')]C',iy), the resulting 
^gi(x) either commutes with the operators to its right, or, in moving to the right of 
those that do not commute with it, produces still other terms that have a gQ\{r) on the 
extreme right-hand side. We therefore observe that the most general expression for the 
commutator of Hq\ and ^q(^) 

Hou e^aiW] = / rfrx"^(r,x)6;a,(r) (80) 

where x°''^{^j ^) is a nonlocal functional of x and r. Although Qq^ (x) is not time-independent, 
it is nevertheless true that a state |^) for which Qq\{^) = for all values of x and a 
will time-evolve so that ^q(x) exp{—iHQ\ = as well. 



The Hamiltonian given in Eqs. (|50D -(p9D relates QCD in the temporal-gauge to its 
Coulomb-gauge formulation, ^^gi is the Hamiltonian of QCD in the temporal gauge. It 
is expressed in terms of gauge-invariant fields, and is displayed in these equations in a 
representation obtained in Ref.[^ by unitarily transforming the standard form of that 
Hamiltonian. But that does not change the fact that Hq\ is still the temporal-gauge 
Hamiltonian. In principle, Hq\ could be used to derive the equations of motion of both, 
the gauge-invariant and the standard gauge-dependent temporal gauge fields, although 
the calculation leading to the latter would be clumsy and exceedingly tedious in this 
representation. In such calculations, we could not neglect the contributions made by 
commutators of operator- valued fields with Hg, since the latter is an indispensable part 
of the temporal-gauge Hamiltonian in this particular representation. In reference to Eqs. 
(|50|)-(|59|), Hq\ and Hg — the two constituent parts of ^^gi — can be given the following 
interpretation: i^Gi can be recognized as a representation of the QCD Hamiltonian in 
the Coulomb gauge, similar to the Coulomb-gauge Hamiltonians obtained by methods 
different from ours. |^, ^, P] These works differ among themselves and from ours 

in various ways — in operator order, in the inclusion or omission of the longitudinal 
component of the chromoelectric field, in whether the chromoelectric field and the glue 
color-charge density are hermitian — but the Hamiltonians given by these authors are 
markedly similar to our Hq\. Furthermore, Aq^^ and Hq^ respectively have the same 
commutation rules with each other, and the components of Hq^- among themselves, as do 
the corresponding Coulomb-gauge fields in Ref. [^], modulo operator ordering. 

Hg has no role in the time evolution of physical systems, since it annihilates all states 
that implement Gauss's law (defined by QQ\{r) |^) = 0). When annihilates a state 
and thus identifies it as implementing Gauss's law, it also annihilates exp{—iHQ\t)\'^) . As 
we have shown, if Hg appears within a sequence of operators constructed from component 
parts of ^^Gh such as might be found in an expansion of exp{—iHQ\t) or in a transition 
amplitude, and if this operator sequence acts on a state that satisfies ^q(x)|\1/) = 0, 
then the resulting expression must vanish, can be commuted through the other oper- 
ator to its right, producing further operators, each of which has a on its right-hand 
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side, until, finally all resulting terms have a Qq^ operator for some value of a and some 
spatial argument, acting on with the result that the original expression {i.e. the 
sequence of operators acting on vanishes. Finally, we note that Hq\ is not affected 
by the operator-ordering problem associated with the direct quantization of QCD in the 
Coulomb gauge referred to in Ref. 0. The gauge- invariant field and negative chromoelec- 
tric field, Aq^^ and Hq^ respectively, obey commutation rules derived from their structure 
in terms of standard gauge-dependent temporal-gauge fields, whose commutation rules 
are simple enough to make the commutation relations of and Hq^, as well as their 
positions in the QCD Hamiltonian, well-defined; 

A remarkably similar state of affairs obtains in QED. When QED is formulated in the 
temporal gauge, and a unitary transformation is carried out which is the Abelian analog 
to the one that takes QCD in the temporal gauge from the C to the Af transformation as 



discussed in Ref. the following result is obtained:[26, 27] The QED Hamiltonian in the 



temporal gauge, unitarily transformed by the Dirac transformation, pm can be described 



as 



Hqed = Jdr [in,(r)n,(r) + ii^^ (r)F,,(r) + ^t(r) {(3m - taA) ^(r) 

- fdrAP{r)Mr) + /rfrrfr'f^^ + H, (81) 
J J 87r|r — r I 

(T) 

where Al designates the transverse Abelian gauge field, and Hg can be expressed as 

Hg = -\jdv (^a,n,(r)^jo(r) + Jo(r)^5,n,(r)) , (82) 

so that all the operator- valued fields that appear in Hqed are gauge invariant, and so that 
Hqed also consists of two parts: the Hamiltonian for QED in the Coulomb gauge and Hg, 
which has no affect on the time evolution of states that implement Gauss's law (which, in 
analogy to the non- Abelian ^Q(r) |\E') = 0, takes the form (9jHj(r)|$) = after the Dirac 
transformation has been carried out).0 When the time derivative is defined as iIHqed 
the equations of motion of temporal-gauge QED result. But since Hg can have no effect 
on the time evolution of state vectors that implement Gauss's law, the only part of Hqed 
that time-evolves such state vectors is the Coulomb-gauge Hamiltonian. This remarkable 
parallelism between QCD and QED prevails because Vc, which we constructed in Ref.[|l|, 
is the non-Abelian generalization of the transformation which Dirac constructed for QED. 



®In Rcfs. 1 26, [2^, 9^11^ is separated into positive and negative frequency parts, so that the gauge- 
invariant charged states is manifestly normahzable. The unitary transformation in these references there- 
fore differs somewhat from the one introduced by Dirac, but has essentially the same effect. 
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3 Multiple solutions of the resolvent field equations 
and Gribov copies 

The fact that we can identify the Coulomb-gauge field with the gauge-invariant fields we 
have constructed, [Q enables us to use previously reported results about the topology of 
the gauge-invariant gauge field to observe how the Gribov ambiguity manifests itself 
in the temporal-gauge formulation we have been discussing. 

The Gribov ambiguity — the existence of multiple copies of fields that obey the 
Coulomb gauge condition — complicates the quantization of QCD in the Coulomb gauge. 
This complication manifests itself as an inability to uniquely invert the Faddeev-Popov op- 
erator Didi. Gauss's law in the Coulomb gauge can be written in the form DidiA^ = —j^, 
in which the Faddeev-Popov operator takes account of the fact that the "total" color 
charge- density (including quark and glue color contributions) is Jg = ge°'''''A'^Il^ + j^. It is 
necessary to invert DidiA^ = — Jq in order to replace Aq in the Coulomb-gauge Hamilto- 
nian with its equivalent in terms of unconstrained field variables. This inversion produces 
not only operator-ordering problems; it also leads to the Gribov problem. Inverting 
DidiAQ = — Jq formally produces the equation Aq = —{Didi)~^jQ, the non-Abelian analog 
of Aq = — (9~^jo in QED. However, whereas is trivially invertible for suitably chosen 
boundary conditions, and manifestly commutes with the Abelian gauge field as well as 
with its conjugate momentum, the differential operator D^di is not similarly invertible 
and does not commute with the non-Abelian gauge field or with the chromoelectric field. 
Aq therefore cannot be uniquely replaced by an expression that depends on unconstrained 
field variables only. 

Various authors have responded to this complication in different ways. Some have 
concluded that the difficulties associated with the Coulomb gauge make it preferable to 
use an axial gauge, such as the temporal gauge, in which, they state, the Gribov problem 



does not exist. Other authors have taken the point of view that, since it is likely 
that Gribov copies of the original gauge field belong to different topological sectors, one 
can ignore the Gribov copies by limiting oneself to the perturbative regime. It has 



also been questioned whether the Gribov ambiguity is truly absent when axial gauges are 
imposed, or whether it always affects non-Abelian gauge theories, but is more obscure 
and harder to detect in the case of axial gauges. Singer has proven that when the path 
integral of a non-Abelian gauge theory is defined over a Euclidean 4-dimensional sphere, 
ambiguities appear regardless of the gauge condition; but, as Singer has noted, in the 
absence of such boundary conditions, there are gauges — in particular axial gauges — to 



which his proof does not apply. |5D|] Others have argued that given these circumstances, it 



is still uncertain whether Gribov ambiguities are endemic to non-Abelian gauge theories 



or whether they only mark the Coulomb gauge as an unfortunate choice of gauge. [P9 



In previous work|I|, |^, ^ Q and in previous sections of this paper, we have given 
series representations of the resolvent field as well as of the inverse of the Faddeev-Popov 
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operator and of a number of Dirac commutators. We have also shown that even when 
functional of gauge fields are given a series representation, that does not necessarily 
eliminate the possibility of representing them nonperturbatively. These functionals can 
be related by integral equations, and nonperturbative solutions of these equations can 
be obtained. In this way, the origin of multiple solutions of gauge fields can be better 
understood 

In Ref.0], we made an ansatz, representing the gauge field in the temporal gauge and 
the resolvent field as functions of spatial variables that are second-rank tensors in the 
combined spatial and SU(2) indices; except in so far as the forms of A]{r) and A]{r) 
reflect this second-rank tensor structure, they are represented as isotropic functions of 
position. In this way, we can represent the longitudinal part of the gauge field in the 
temporal gauge as 



a: 



7 L/ 




and the transverse part as 



TB(r) + e 



i-yn—Tc{r) 

r 



^3) 



(84) 



where N'{r), T^ir), Tsir) and Tc{r) are isotropic functions of r, the prime denotes dif- 
ferentiation with respect to r, and the transversality of AJ'^iv) requires that 



dr 



(TTa 
dr 



0. 



^5) 



An entirely analogous representation of the resolvent field enables us to relate it to the 
gauge field through the nonlinear integral equation described in Section |l]. As a result of 
this analysis, we have been able to show that it is possible to represent the resolvent field 
asi 

rir^\ (Jf \ r„ 
*7 ) I ~ + '^^ ) + ^n^TV^c- (86) 



5r 



where 



and 



gr 



gr 



'Mcos{N + AT) - sin(Ar + M) 



"Arsin(Ar + U) + cos(F + M) 



cos(Ar + AT) - ll - Tc sin(Ar + M) (87) 



+ Tc [cos(Ar + U)-1\+Ta sin(Ar + U). 

(88) 
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Similarly, the gauge-invariant gauge field can be expressed as a functional of Af and of 
Af, Ta, Tb and Tc as shown by 



1 

gr 



-J 7 n 



cos(Ar + A/") - 1 + Arsin(Ar + M) 



+ ( 5.. - ^ 



'Nco^iN + AT) - sin(Ar + M) 



r dr 



+rA 



COS 



+Tc |e,^„^cos(Ar + AT) - (^6,^ - ^) sin(Ar + AT) } . 



^9) 



With these representations, the nonlinear integral equation that relates the resolvent field 
to the gauge field (in the temporal-gauge) is transformed to the nonlinear differential 
equation 



'Arcos(Ar + AT) - sin(Ar + AT)' 



d^Af djf 2 
du^ du 

[Ta [cos(Ar + J\f)-l]-Tc sin(Ar + AT) } = (90) 

where u = ln(r/ro) and tq is an arbitrary constant. AT — the dependent variable in Eq. 
(|90|) — can be directly related to the resolvent field A] by 



(91) 



Similarly, can be related to the gauge field in the temporal (Weyl) gauge as shown by 

^ = ^7|^"W- (92) 
~|) and by the transverse part of the gauge field 



Ta, Tb and Tc are determined by Eq. (| 
as shown by 

Ta + Tb and by ^e'^^^^A] ^ = Tc . 
2 r 



(93) 



In solving Eq. (PO]), boundary conditions are imposed on AT and on the source terms Af, 
Ta, Tb and Tc- Ta,Tb and Tc are required to vanish as m ^ ±oo; A/" is required to vanish 
as M— s> — oo, and to either vanish or approach specified limits as u-^oo. Af is required to 
be bounded in the entire interval — oo<m<oo. For convenience, we normalize AT so that 
AT^O when — oo. 

One of the questions addressed in Ref.[Q is the following: given a specified set of values 
Af{r), Ta^v), Tb^v) and Tcir), what behavior is possible for AT, if AT is required to be 
bounded in the entire interval 0<r < oo and to obey Eq. (pO])? In posing this question. 
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77 is seen to be the unknown function, and Af{r), T^ir), Tsir) and Tc{r) are source 
terms that drive the solution, subject to the required boundary conditions. In Ref. 0, 
we presented numerical integrations of Eq. (^) which demonstrated that a number of 
different solutions of this equation — and, hence, a number of different — can be 
based on the same gauge-dependent gauge field characterized by a single set of values 
of A/", T4, Tb and 7^. A related question — what form for Ag, -(r) corresponds to the 
different possible functional forms of A/" — can also be posed. When the gauge field 
Aj = 0, the resolvent field A] and the gauge-invariant field need not vanish. In that 
case, Eq. ( pOD reduces to the autonomous^ "Gribov equation" |l^ 



-r^ + 2sm(Ar) = (94) 

which is also the equation for a damped pendulum with AT representing the angle with 
respect to the pendulum's position of unstable equilibrium, and u representing the time.|^ 
In the application of this equation to the Gribov problem, AT must remain bounded not 
only in the interval 0<m < 00, but also in the larger interval —00 < n < 00 to include the 
entire configuration space 0<r < 00. 

The restriction that AT must be bounded in the entire interval (— oo<m<oo) severely 
limits the allowed solutions of Eq. (p^). For solutions of Eq. (p^) that are bounded in the 
entire interval (— oo<m<oo), there is a single, unique, phase plot. One branch extends 
from the unstable saddle point (corresponding to — cx)), at which AT = 0, to a stable 
point (corresponding to m^oo), at which A/" = vr. The other branch extends from the 
same saddle point (at which u corresponds to —00) to a stable point at AT = — tt, also 
corresponding to n^oo. The two branches are identical, except that for each point on one 
branch the values of AT and dH jdu correspond to values —AT and —dJJjdu on the other. 
Numerical solutions are obtained by setting AT and dJJ j du equal to the same small value 
at some large negative value of m, in order to discriminate against solutions that become 
unbounded as u approaches the saddle point as — 00, as discussed in Ref.|Q. We can 
choose different large negative values of u at which to set AT = dJJ/du = 0; if, in one case, 
we choose Ua and in another Ub, we obtain solutions Uaiu) and J\fb{u) respectively, where 
J\fb{u) = J\fa{u+\J) and U = Ub — Ua- Similarly, changing the magnitude of the small value 
of Uiua) = dNiua) / du = £ to J\f{ua) = dN{ua)/du = e' also has the effect of shifting 
the functional form of Af{u) from u to u + Uq for a particular value of Uq. Replacing the 
initial condition J\f{ua) = dj\f{ua) /du = e, for a small e, with J\f{ua) = dJ\f{ua) /du = — e, 
has no other effect than changing the signs of A/'(m) and dSf /du. Figure 1 illustrates some 
of these relationships. 

To study the variation in the form of AT that is allowed when AT is represented in 
configuration space in the form A/'(M(r)) = A7'(ln(r/ro)), we observe the following: For the 

^Eq. ( p4| ) is autonomous because the variable u does not appear explicitly in this equation; it appears 
implicitly only as the argument of A/" and of its derivatives. 
^°the function a in Ref.p^] is related to AT by AT = 2a. 
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shift J\fb{u) = Uaiu + U) with U represented as U = ln(ro/-Ro) for an appropriate i?o, 

ATfe (ln(r/ro)) = Jfa (ln(r/ro) + ln(ro/i?o)) = Jfa (ln(r/i?o)) . (95) 

It is clear that when AT is required to be bounded in the entire interval (0<r<oo), a 
change in the constant tq and an overall change of sign are the only changes allowed in 
AT (ln(r/ro)). In particular, in every such case, Ar(r = 0) = and Ar(r— >-cxd) = ivr. When 
the restrictions on A/", T4, Tb and Tc that led to Eq. (p^ are apphed to Eq. (^), we 
find that 

Ku(r)](o) = ^^..n^, (96) 

where the subscript (0) designates the solution that corresponds to the "pure gauge" case 
in which the gauge-dependent gauge field AJ vanishes. [^Gii('')]{o) can easily be recognized 
as a "hedgehog" solution. 



Some authors have suggested that Eq. (|94D — the Gribov equation — has a variety 
of solutions for which A7'(r— i>oo) can be any integer multiple of vr, [|32| and that differ- 



ent integer multiples correspond to different topological sectors connected by large gauge 
transformations. This suggestion, which is motivated by the model of a damped pen- 
dulum, neglects the fact that in the Gribov equation u must be bounded in the entire 
interval {—oo<u<oo), and that the only solutions of the damped pendulum problem that 
can remain bounded in this entire interval as the time u is extrapolated backwards to —00, 
are those for which the pendulum initially is at rest in its unstable equilibrium position. 
And, with this initial position, the damped pendulum is unable to execute multiple turns 
before coming to rest at equilibrium in a stable configuration. 

Gribov explicitly noted that the necessity of requiring the solutions to Eq. (|94D to be 
bounded in the entire interval —00 < m < 00 limits the asymptotic values of A7'(m— >oo) 
to ±7r.[Q He therefore also considered the transverse gauge field that results when an 
arbitrarily chosen transverse gauge field is gauge-transformed, so that the new field is not 
pure gauge, but is given by 

= UkiU-^ + lUd.U-^ (97) 

-i(j)(r) and Aj = A"^— with /K'r = e'^^^fh). 

2r J 2 r'^ 

A^ and Aj both are transverse and therefore belong to the Coulomb gauge. The transver- 
sality of A^ leads to the equation 

^ + ^_2,,„(^),l_;(.)).0 (98) 

which is not autonomous, and does have the multiple solutions that correspond to Gribov 
copies belonging to different topological sectors. 
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In our work, Eq. (^) — of which Eq. (0) is a special case — also has inhomogeneous 
source terms that prevent it from being autonomous. But, unlike Eq. (]98|), Eq. (0) is 
not obtained by gauge-transforming from one Coulomb-gauge field to another. It is a con- 
sequence of the transformation from "standard" gauge-dependent fields in the temporal 
gauge to the corresponding gauge-invariant fields. Our Eq. (|9^) describes this trans- 
formation for the case that the initial gauge-dependent temporal-gauge field is set = 0, 
which clearly corresponds to the "pure gauge" case — Gribov's Eq. ( pSD with f{u) = 0. 
But, more generally, the relation between our Eqs. (|90D and (0) is different from the 
relation of Gribov's Eq. (|98|) to that same equation with f{u) = 0. Our Eq. (^) would 
not reduce to an equation of the form of Eq. ( ]98| ) even if the initial gauge- dependent field 
were chosen to be purely transverse. Eq. (pOD shows that, in that case, there would be 
the additional "source"-term 2(yfro exp(M)7^(cosAr— 1) (where AT corresponds to in Eq. 
(|98|) ). Nevertheless, in spite of these differences, the fact that we are able to identify the 
gauge-invariant gauge fields we constructed with the Coulomb-gauge fields enables us to 
interpret the multiple solutions of Eq. (|90D, which we demonstrated in Ref. 0, as Gribov 
copies. These multiple solutions represent Gribov copies of the gauge- invariant gauge 
fields in the temporal gauge, which manifest themselves even when there are no such 
copies, and no ambiguity, in the case of the "standard" gauge-dependent temporal-gauge 
fields. 



4 Discussion 

In this work, we have established a relationship between the gauge-invariant temporal- 
gauge and the Coulomb-gauge formulations of two-color QCD: The temporal-gauge QCD 
Hamiltonian, when represented entirely in terms of gauge-invariant operator-valued fields, 
is not identical to the Coulomb-gauge Hamiltonian represented by Hq\ given in Eq. (W^). 
But the two are physically equivalent — i. e. they lead to identical values of observable 
quantities within the space of states in which Gauss's law has been implemented. More- 
over, we have shown that the gauge-invariant temporal-gauge fields obey commutation 
rules that are the same as those for Coulomb-gauge fields, modulo operator-ordering ambi- 
guities, when the longitudinal components of the Coulomb-gauge chromoelectric field are 
retained, as they are in our work and in Schwinger's treatment of the Coulomb gauge. ||^ 

We have also shown, through reference to specific numerical calculations, that there 
are Gribov copies of gauge-invariant fields in the temporal-gauge formulation of QCD, 
even though there are no Gribov copies of the gauge-dependent temporal-gauge fields. We 
have demonstrated that the gauge-invariant fields in the temporal gauge obey a nonlinear 
integral equation, which — subject to an ansatz — can be transformed to a nonlinear dif- 
ferential equation that has multiple solutions corresponding to a single gauge-dependent 
field. These solutions must be bounded in the entire configuration space. Thus, the non- 
linear differential equation that embodies the imposition of Gauss's law and the implemen- 
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tation of gauge invariance in the temporal gauge, and the requirement of boundedness, 
lead to the multiple solutions that can be identified as Gribov copies. This is consistent 
with our demonstration of a close resemblance between the gauge-invariant formulation 
of the temporal gauge and the Coulomb-gauge formulation of QCD. It is also consistent 
with the fact that the nonlinear differential equation that relates the gauge-invariant and 
the gauge- dependent temporal-gauge fields — Eq. (pO]) — has a form very similar to the 
one that Gribov used to demonstrate the non-uniqueness of Coulomb-gauge fields — Eq. 
(|98|). 



We are therefore led to conclude that the reason why the Gribov ambiguity does not 
arise when QCD is quantized in the temporal gauge — and most authors who discuss the 
quantization of QCD in the temporal gauge either do not mention the Gribov ambiguity 
or state that there is no Gribov ambiguity in the temporal gauge |^ — is that, in 



sharp contrast to the Coulomb gauge, quantization in the temporal gauge proceeds to 
completion without requiring the imposition of Gauss's law. The Gribov copies arise in 
the temporal-gauge formulation only after the theory has been quantized, and Gauss's law 
then is implemented and gauge-invariant fields are constructed. Our results are consistent 
with the conclusion that the Gribov ambiguity is a fundamental attribute of non-Abelian 
gauge theories. Gribov copies do not arise when a non-Abelian gauge theory is originally 
quantized in the temporal or other axial gauges, because, then. Gauss's law remains 
unimplemented. The Gribov ambiguity does manifest itself in these gauges, but only 
with the introduction of the gauge-invariant fields and the imposition of Gauss's law. 
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Appendix 

In this Appendix, we will evaluate the quantities Vui and TZyb defined in Eq. and 
show that they have the required properties. 

For the case of Vui, we evaluate VcdiV^^, where Vc is most conveniently represented 
as shown in Eq. (0), so that 



Xi 



Vcd.Vc' = exp (-^(^Z"(r)^) d.exp (^(?Z"(r)^) . (99) 



For = gZ"^ and $ = we obtai 



am 



Xi = ^ cos ( ^) - i^-^ sin (^1^ I di jcos (^^^ + i^-^ sin (^^]\. (100) 
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Eq. then determines that 

. + (£|£ _ (i - + ) . (101) 



Similarly, Eq. (p4D also determines that 7?.fea is given by 



7^,, = 5,, cos $ + <l>'^<l>'' j + e"''^^'^ ^-^j . (102) 

Inspection of Eq. ( |102| ) demonstrates that simultaneous substitution of — $ for $ and 
exchange of the subscripts a and b in 7?.f,a leaves the expression on the right-hand side of 
this equation unchanged, thus proving Eq. (|39D . 
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Figure 1: 



Figure Caption 



Three solutions of Eq. (Q) with AT bounded in the interval (— oo<m<oo). Of the two 
solutions with Ar>0, one is obtained with AT = dN /du = 10^^ at u = —11 and the other 
with AT = dJI /du = 10~^^ at M = —11. The plot for the latter solution is identical to the 
plot of the former shifted to the right on the u axis. Both solutions approach TV" = vr as 
u^oo. The solution with J\f<0 is obtained with AT = dj\f/du = — 10~^ at m = —18. In 
this solution, A/"— — vr as m— >oo. The three solutions are precisely identical modulo shifts 
along the u axis and reflection in it. 
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